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Abstract. The contribution of multiple forward scattering in Coulomb focusing of 
low-energy photoelectrons at above-threshold ionization in mid-infrared laser fields 
is investigated. It is shown that the high-order forward scattering can have a 
nonperturbative effect in Coulomb focusing. The effective number of rescattering 
events is defined and is shown to depend weakly on laser intensity and wavelength. 
Nevertheless, the relative contribution of forward scattering in Coulomb focusing and 
the Coulomb focusing in total decrease with increasing laser intensity and wavelength. 

1. Introduction 

The Coulomb field of the atomic core can play a significant role in the strong-field 
photoionization process essentially modifying the dynamics of low-energy electrons. It 
is responsible, in particular, for the appearance of a rich structure in the momentum 
distribution of photoelectrons near the ionization threshold [IIEIEIIIIEJEJITIIE], for 
frustrating the tunneling ionization jH] and for the creation of a low-energy structure 
(LES) in photoelectron spectra in mid-infrared laser fields [TO], [HI [121 021 EH]- The 
Coulomb field focuses low-energy electrons towards the laser polarization direction 
which is mostly due to multiple rescattering [JS] of ionized electrons by the atomic 
core at large impact parameters and is termed Coulomb focusing (CF) [161 [TTl [18]. 
For a theoretical description of Coulomb field effects, different modifications of the 
strong field approximation [19] have been developed [U EOj EU [221 E3]. In addition, 
the classical trajectory Monte Carlo (CTMC) method has been successfully employed, 
see e.g. [211 ESI HH1 ttH EE], for estimation of the effects which are not intrinsically 
quantum mechanical. 

Recently, the strong field physics in mid-infrared laser fields has attracted a lot of 
attention in connection with the possibility of improving high-order harmonic generation 
with mid-infrared driver fields [2b] . In mid-infrared laser fields, when the Keldysh 
parameter is small 7 = \Jl p /2U p <C 1, the electron dynamics after tunneling is mainly 
classical. This is because the characteristic energies of the process, I p and U p , greatly 
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exceed the photon energy in this regime u <C I p <C U p . Here, I p is the ionization 
potential, U p = Eq/Au 2 the ponderomotive energy, E and u are the laser field amplitude 
and frequency, respectively (atomic units are used throughout). In this regime, the 
classical features of the three-step model [15] are conspicuous and not obscured by 
interference effects. Two recent experiments by Blaga et al. [TU] and Quan et al. [TT] 
on the photoionization of atoms and molecules in strong mid-infrared laser fields reveal 
a characteristic spike-like LES in the energy distribution of electrons emitted along the 
laser polarization direction. The CF is responsible for the effect p2J E2 UJ- More 
concretely, the LES arises due to multiple forward scattering (FS) by Coulomb field 
rT5] . The CF is usually predicted to decrease with increase in the laser intensity and 
wavelength because the average rescattering velocity and the impact parameter increase 
in such circumstances. As a consequence, one may expect that the contribution of 
high-order FS should also decrease. In this context, it was surprising that at a large 
wavelength of the mid-infrared laser field, the multiple FS plays a decisive role for the 
creation of the LES. 

In this paper, we investigate how the contribution of different components of CF 
depends on laser intensity and wavelength. Our investigation is limited to the classical 
interaction regime in mid-infrared laser fields. Separate components of CF are identified 
which scale differently with laser parameters: CF which happens immediately after 
ionization - initial CF (ICF); CF due to the electron FS on recollision with atomic core, 
and asymptotic CF (ACF) when the electron momentum is disturbed by the Coulomb 
field after the laser pulse is switched off. Special attention is devoted to the contribution 
of the high-order FS events and to the definition of the effective number of FS events. 
We use the CTMC method with tunneling and the Coulomb field of the atomic core 
fully taken into account. 

2. The method 

The 3D CTMC method employed in this paper is developed as follows. (1) An ensemble 
of electrons is formed corresponding to the tunneled electron wave packet according to 
the Ammosov-Delone-Krainov (ADK) theory [27]. The electrons are born with the 
following initial conditions. The electron initial position along the laser polarization 
direction is derived from the Landau's effective potential theory [28]. The transverse 
coordinates of the initial position are zero. The initial longitudinal momentum is zero 
and the transverse one follows the corresponding ADK distribution [29|. (2) The 
electron wave-packet propagates in the field of a laser pulse and Coulomb potential 
via the solution of Newton equations. (3) The positions and momenta of electrons 
when the laser pulse is switched off are used to calculate the asymptotic momenta at 
the detector [30]. (4) Each trajectory is weighted by the ADK ionization rate and the 
initial transverse momentum distribution function [29]. (5) The shape of the laser pulse 
is half-trapezoidal: For the first ten cycles, the field has a constant amplitude and is 
ramped off within the last three cycles. The electrons are launched within the first half 
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cycle since there are no multi-cycle interference effects in the classical calculation. The 
ensemble consists of 10 6 particles and the convergence is checked via double increase of 
this number. The target atom is neon with ionization potential I p = 21.56 eV which 
can endure a maximum laser intensity Iq 8.66 x 10 14 W/cm 2 |31j. The process is in 
the tunneling regime, e.g., 7 ~ 0.2 at Iq = 7.24 x 10 14 W/cm 2 and wavelength A = 2/im. 

We estimate the contributions of the multiple FS, ICF and ACF to the total CF 
in the following way: (1) For each rescattering event at the moment t s , the minimal 
distance from the core r s , the distance from the core in the transverse plane (with respect 
to the laser polarization direction) p s and the electron momentum p s are determined 
numerically. Then, the transverse momentum change Sp± due to the Coulomb potential 
V(r) at the s-th forward scattering event is estimated as 5p± s ~ / V ±V(r(t))dt ~ 
— (p s /r^)5t s , where 5t s is the rescattering time duration. When the electron velocity 
p s in the FS event is large, 5t s ~ 2r s /p s . In the opposite case, 5t s ~ 2^2r s /\E(t s )\ is 
determined by the laser field E(t s ) at the s-th rescattering moment t s . Accordingly, 

5p ±s = - 2p s /(r 2 p s ), if p 2 s > r s \E{t s )\ 

5p ±s = -2 3 ^ 2 p s /^ri\E(t s )l otherwise. (1) 

(2) The transverse momentum change due to ICF is estimated numerically as the 
deviation of the exact transverse momentum from that neglecting the Coulomb 
potential, after a half laser period following the ionization moment tf. dp^P = p±(ti + 
T/2) — p± (t{ + T/2), where T is the laser period and p±* is the electron transverse 
momentum neglecting the Coulomb field. The numerical estimate for 8pf is slightly 
larger in absolute value than the analytical one [7]: 

5p ( [ ] ^-2p l± \E(t t )\/(2I p ) 2 , (2) 

with the initial transverse momentum pij_. (3) We estimate the ACF contribution via 
numerical comparison of the asymptotic electron momentum with the one after switching 
off the laser pulse. 

3. The results 

The CF is mainly due to multiple small-angle scattering. It is significant only for 
low energy photoelectrons, and we will examine the dynamics for such electrons in 
details. The CF is characterized by the transverse momentum change Sp± induced by the 
Coulomb field [13] which depends on the ionization phase ipi = uU. We restrict ourselves 
to ionization phases and to trajectories which contribute to the low energy part (up to 40 
eV) of above-threshold ionization spectrum emitted along the laser polarization direction 
within an opening angle of ±2.5°. The electrons, which are emitted out of the laser 
polarization direction, have experienced large-angle scattering, their CF is interrupted 
and, consequently, their dynamics is not typical for CF. The laser intensity dependence 
of the different CF components is shown in Fig. [T] and the wavelength dependence 
in Fig. [2] For each (pi, the transverse momentum change is shown for the electron 




Figure 1. (color online) The transverse momentum change 8p± versus the ionization 
phase ifi — ioti. The CTMC simulation for a neon atom in a mid-infrared laser field 
with a wavelength A = 2/im for the following laser intensities: (al-a3) I = 1.81 x 10 14 
W/cm 2 , (bl-b3) / = 3.62 xlO 14 W/cm 2 and (cl-c3) / = 7.24 xlO 14 W/cm 2 . (al,bl,cl) 
The total transverse momentum change (marked as "exact" ) , the estimation of ICF and 
ACF as well as of the total transverse momentum change (marked as "ICF+ ACF+FS" ) 
as described in Sec{2] (a2,b2,c2) Sp± due to the s-th order FS events (s is indicated in 
the inset) and due to ICF. (a3,b3,c3) The ratio of the 6p± at the s-th order FS events 
to the first-order one. The ionization phase (p^' corresponding to the threshold of the 
multiple FS is marked by an arrow in (al). The peak and the plateau of 5p± for the 
2nd FS are marked by arrows in (a2). The maximum of the laser field is at ifi = Tt/2. 



trajectory which has the maximal probability among the contributing trajectories at 
this ionization phase. We calculate the total transverse momentum change exactly via 
the CTMC simulation, see the curves marked as "exact" in Figs, [I] (al,bl,cl) and [2] 
(al,bl,cl). Further in Figs. [l]and[2j we show the results of the estimate of Sp± due to 
the s-th FS (s < 6), ICF and ACF as described in Secjl} To show the accuracy of our 
estimations, we sum up all contributions to 5p± and comapare it with the exact result, 
see the curves marked as "ICF+ACF+FS" in Figs. [I] (al,bl,cl) andg (al,bl,cl). 

From the analysis of Figs. [T] and [2j the following conclusions can be drawn. First 
of all, Figs. [T] (al,bl,cl) and [2] (al,bl,cl) show that the contribution of ACF to the 
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Figure 2. (color online) The transverse momentum change 8p± versus the ionization 
phase ifi . The CTMC simulation for a neon atom in a mid-infrared laser field with an 
intensity of I = 1.81 x 10 14 TT/cto 2 for the following wavelengths: (al-a3) A = 2 /im, (bl- 
b3) A = 2.828 /im and (cl-c3) A = 4/j,m. (al,bl,cl) The total transverse momentum 
change (marked as "exact"), the estimation of ICF and ACF as well as of the total 
transverse momentum change (marked as "ICF+ACF+FS") as described in Sec{2] 
(a2,b2,c2) Sp± due to the s-th order FS events (s is indicated in the inset) and due to 
ICF. (a3,b3,c3) The ratio of the 5p± at the s-th order FS events to the first-order one. 



total 5p± is generally negligible (the main contribution is at ionization phases near the 
peak of the laser field within the ionization phase interval of 5ipi « 0.02). It decreases 
with increasing intensity and does not change with wavelength. The contribution of 
ICF to the total 5p± increases with increasing intensity and remains almost constant 
with wavelength which is consistent with the estimate of Eq. ([2]). The contribution 
of ICF still constitutes a small fraction of the total Sp± (less than 10%) for ionization 
phases 7r/2 < ifi < ipf* (the maximum of the laser field is at ipi = vr/2), where multiple 
scattering takes place but competes with the single scattering contribution at ifi > </?j , 
especially at high intensities and wavelengths (see the estimate below, Eq. &). The 
ionization phase (p\ marks the threshold of the multiple FS, see the indication of (p[ 
in Fig. bl (al); (p^ ~ 1.8 at a laser intensity of 1.8 x 10 14 W/cm 2 but decreases slightly 
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Figure 3. The electron trajectories at ipi = 1.75 (thick line) and <pi = 1.64 (thin 
line). The x-coordinate along the laser polarization direction is scaled by Eq/cj 2 and 
the transverse z-coordinate by pij_/uj. The position A corresponds to the plateau and 
B to the peak of the transverse momentum change at the second forward scattering. 



with increasing intensity. 

The main CF contribution to the total 5p± comes from FS which determines the 
shape of the curve 5p± versus the ionization phase, shown in Figs, [I] (al,bl,cl) and 
[2] (al,bl,cl). The transverse momentum change 8p± s due to the s-th order FS in the 
case of different laser intensities and wavelengths are shown in Figs. [I] (a2,b2,c2) and [2] 
(a2,b2,c2), respectively. 5p± s has a characteristic dependence on the electron ionization 
phase which is qualitatively the same for each scattering order. The 5p± s increases 
sharply with decrease in the ionization phase from the threshold value (it is different for 
different scattering orders), reaches the peak and then decreases slowly down to a flat 
plateau, the latter having an increasing tail with the further decrease of the ionization 
phase. Although the contribution of FS decreases on average with increasing order, a 
higher-order FS can make a larger contribution in some phase intervals than a lower- 
order one. 

Let us estimate the values of the peaks and plateaux of the transverse momentum 
change due to high-order FS . The peak in the Sp± s for the s-th order FS arises when the 
electron trajectory touches the z-axis at a recollision (the coordinate center is chosen at 
the atomic center, x-axis is in the laser polarization direction and z-axis in the transverse 
direction, see Fig. [3]). In this case, r s w p s , E(t s ) ~ E , the electron momentum p s is 
nearly zero, see Fig. [5] (a-d), and 5p± s is determined by the second expression of Eq. 
([!]). The ionization phases corresponding to the peak and the plateau of the s-th order 
FS are indicated in Fig. [4] by (pfj and <pfs \ respectively. The impact parameter at the 
peak of the s-th order FS can be estimated as pP eak ~ Pi±7r(s + l)/ui, see Fig. |3j Then, 
using Eq.Q we have 



The plateau in the 5p± s corresponds to the FS case when the electron velocity in the 
FS point is the largest, see Fig. p] (c). The latter is determined by the amplitude of 




(3) 
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Figure 4. (color online) The parameters of the s-th order FS (s is indicated in the 
inset): (a) p s the distance from the atomic core in the transverse plane, (b) r s the 
distance from the core, (c) p s the momentum, (d) E(t s ) the field at the scattering 
moment. ip\ p J and ip^P are the ionization phases corresponding to the peak and 
plateau for the s-th order FS, respectively, which are indicated with dashed lines. The 
CTMC simulation for a neon atom in a laser field with I = 1.81 x IQ^W/cm 2 and 
A = 2 /mi. 



the velocity oscillation in the laser field: v s ~ (3 s Eo/u, with (3 S « 0.8 at even s and 
p s ~ 1 at odd s, according to the numerical results. The impact parameter in this case 
is estimated as pP lateau ~ pi^irs/u, which yields 

«| (4) 

The estimates of Eqs. ([3]) and Q, which are in agreement with the numerical 
calculations presented in Figs. [T] and |2j show that the peaks for the higher order FS 
((s + l)-th order) can exceed the plateaux of the lower-order FS (s-th order). In fact, 
this ratio is 



which is between 1.1 and 1.2 for s = 2 — 6 at 1.81 x 10 14 W/cm 2 . Especially, the peaks 
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Figure 5. (color online) The CTMC simulation for a neon atom in a mid-infrared 
laser field: (a,c) Photoelectron spectra. (b,d) The transverse momentum change 
versus ifi. (a,b) A — 2/ito for different laser intensities indicated in the inset. (c,d) 
/ = 1.81 x 10 14 W/cm 2 for different wavelengths indicated in the inset. 



of the even order FS (2nd, 4th,...) are larger than the corresponding odd FS plateaux 
(1st, 3rd,...). The plateaux of the even order (2nd, 4th,...) FS are comparable with that 
of the corresponding odd FS (1st, 3rd,...). These are because the velocity at even FS 
events is smaller than that at odd FS. Eq. (|5| shows a remarkable feature that the peak- 
to-plateau ratio increases with increasing intensity and wavelength. In particular, due 
to a larger contribution of the 6th order FS peak with respect to the plateau of the 5th 
order FS at intensity 7.24 x 10 14 W/cm 2 , see Fig. [I] (c2), an additional oscillation in the 
5p± dependence on the ionization phase arises, at (fi « 1.64 in Figj5] (b), which induces 
an additional lower energy peak (at about 9 eV) in the energy distribution within the 
LES, see Fig. [5] (a). 

Using Eqs. ^ and Q, we calculate the ratio of ICF to the plateau of the first-order 

FS: 

^plateau * [ 2 I p U J ' U 

The latter confirms the above statement that ICF can compete with the single scattering 
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contribution to CF at high intensities and wavelengths. For instance, at typical 
parameters Pi± = 0.1 a.u., I p = 0.79 a.u., I = 7 x 10 14 W/cm 2 and A = 2fim the 
ratio in Eq. (6) amounts to Spf 1 /5p± f oau ~ 0.4. 

The contribution of the FS (as well as the total CF effect) decreases with increasing 
intensity and wavelength because of the increased scattering velocity and the impact 
parameter, see Eqs. ^ and Q. For the same reason the curves in the phase space 
move down with higher intensities and wavelengths in Figs. [5] (b) and [5] (d). However, 
the ratio of 5p± of different scattering orders does not significantly variate with variation 
of laser intensity and wavelength. Thus, the peak of the transverse momentum change 
due to the s-th FS scaled with that of the first FS can be estimated 

It shows that the relative role of the s-th FS even incerases slowly with increasing laser 
intensity and wavelength. This can also be seen from Figs. [T](a3,b3,c3) and[2](a3,b3,c3). 
The mentioned feature can be interpreted as a slow variation of the effective number of 
scattering. 

4. The effective number of scattering 

The effective number of scattering N e ff can be defined employing Eqs. ^ and Q. We 
define N e ff as follows. First of all, we choose a small parameter e C 1 to determine 
the accuracy in which the contributions of the high-order FS to the total 5p± can be 
neglected with respect to that of the first FS Sp±\. The effective number of FS N e ff is 
defined to determine the highest order of FS which makes nonnegligible contribution to 
the total d~p±. Namely, at a given e, 

5pi e ?MT U > e, if s < N eff , 

5p p ± e f/SpfT u <e, if s>N eff . (8) 
The latter definition yields the following expression for the effective number of FS 

^^f^y /3 -i. P) 



nupi 



For the intensity range shown in Fig. [TJ our criterion gives for N e ff a number between 
5 and 6 (for the concreteness e = 0.1 is assumed). The number of scattering according 
to Eq.g increases slowly with increase in the laser intensity and wavelength. 

Generally, the total CF effect decreases with intensity and wavelength because the 
main contribution in CF comes from the multiple forward scattering which, in total, 
decreases as expected, see Fig. [5] (b) and (d), where the curves in the phase space move 
down, generally, with increasing laser intensity and wavelength. The exception is the tail 
at ifi > ipp in Fig. [5] (b) which is due to the competition of the first FS with the initial 
CF as mentioned above. The most important property is that the ratio of the transverse 
momentum changes due to FS of different orders are almost constant with increasing 
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laser intensity and wavelength, see Figs. [T] (a3,b3,c3) and [2] (a3,b3,c3). Therefore, 
the shape of the phase space distribution (dependence of 5p± on the ionization phase) 
remains similar, see Fig. [5] (b) and (d). Note that the consecutive slope changes of the 
phase space distribution are responsible for the creation of LES and are determined by 
the contributions of the second, third and forth FS [32]. 

We point out some features which distinguish the wavelength dependence of CF 
from the intensity dependence, see Fig. [2j The FS contribution to the total Sp± 
decreases more strongly with increasing wavelength as Eqs. ^ and Q indicate. Due 
to the latter, the knee of the phase space distribution becomes less prominent at higher 
wavelengths. The relative contribution of ICF to the total Sp± with respect to FS 
increases with increasing wavelength but to a lesser extent than in the case of the 
intensity dependence. This is because of the Pi± factor in Eq. ^ which decreases with 
increase in the wavelengths. 

5. Conclusion 

At above-threshold ionization in the realm of intensities and wavelengths corresponding 
to the classical regime u <C I p and 7 < 1, multiple forward scattering of an ionized 
electron has a nonperturbative contribution to Coulomb focusing. In some regions of 
ionization phase (photoelectron energy), the contribution of the higher-order forward 
scattering to the total Coulomb focusing can dominate the lower-order one which creates 
local peaks in the photoelectron spectra. The effective number of scattering does not 
depend significantly on laser intensity and wavelength. 
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